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Summary 

The  following  missile  assignment  problem  is  considered*  Missiles 
ore  to  be  assigned  to  targets  in  two  distinct  steps*  First,  each 
mijsile  is  programmed  so  that  it  can  be  fired  at  any  one  of  a  small 
number  of  targets,  the  number  of  targets  being  the  missile  capability. 

The  programming  of  the  missiles  is  represented  by  a  qualification 
matrix  Q.  Second,  if  battle  occurs,  all  missiles  are  to  be  assigned 
to  targets  and  launched.  Each  missile  must  be  assigned  to  a  target  for 
which  it  is  programmed.  It  is  assumed  that  only  a  random  subset  X 
of  the  missiles  will  actually  be  available  for  battle,  and  so  the 
assignment  must  be  made  for  a  reduced  qualification  matrix  Q(X).  The 
questions  considered  are  "what  is  an  optimal  assignment  given  the  reduced 
qualification  matrix  Q(X)?",  and  "what  can  be  expected  from  this  assign¬ 
ment?"  Use  of  a  damage  function  is  proposed.  An  optimal  assignment  is 
one  which  maximizes  the  value  of  the  damage  function.  The  damage  function 


may  be  chosen  to  represent  a  wide  variety  of  optimization  requirements. 

The  main  part  of  the  paper  describes  Monte  Carlo  procedures  for  estimating 
the  expected  damage  and  the  probability  that  the  damage  will  be  at  least 
c  for  any  nimiber  c.  In  an  appendix  several  examples  are  given  to 
illustrate  the  use  of  the  damage  function.  Another  paper  [l]  describes 
an  algorithm  for  finding  an  assignment  which  maximizes  the  damage  function* 
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Introduction 


This  note  deals  with  certain  mathematical  problems  connected 
vd-th  the  assignment  of  missiles  to  targets.  The  problems  are  those 
which  arise  because  the  missiles  have  a  multiple  target  capability. 

The  assignment  of  such  missiles  to  targets  occurs  in  two  distinct 
steps.  First,  each  missile  is  programmed  so  that  it  can  be  fired  at 
any  one  of  a  small  set  of  targets.  The  number  of  targets  for  which 
the  missile  may  be  programmed  is  its  target  capability.  Second,  when 
battle  actually  occurs  and  the  missile  is  to  be  fired,  it  is  assigned 
to  one  of  the  targets  for  which  it  is  programmed.  The  second  step, 
i.e,  the  assignment  of  the  missile  to  a  target,  takes  almost  no  time. 
However,  the  first  step,  the  programming  of  the  missile  for  a  set  of 
targets,  is  very  time  consuming.  Therefore,  we  assiune  that  if  battle 
actually  occurs,  the  missile  cannot  be  reprogrammed  but  must  be  fired 
at  one  of  the  targets  for  which  it  is  already  programmed. 

When  battle  occurs,  it  may  be  that  not  all  of  the  missiles  are 
available.  Still  because  of  the  multiple  target  capability,  there  will 
be  many  different  ways  of  assigning  the  available  missiles  to  targets. 
One  problem  is  to  pick  the  best  possible  assignment.  Another  problem 
is  to  evaluate  the  effectiveness  of  the  original  missile  programming. 
Because  of  the  complexity  of  the  final  assignment  procedure  (the 
procedure  which  picks  out  the  best  assignment  of  available  missiles), 
this  problem  must  be  done  by  the  Monte  Carlo  method. 

To  describe  the  problems  more  precisely  emd  introduce  some  of  our 
terminology,  let  us  look  at  a  simple  example. 
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Consider  a  group  of  nine  missiles,  each  missile  having  a  capability 
of  two  targets,  and  suppose  that  there  are  six  teurgets.  The  group  might 
be  set  up  with  the  missiles  prograinned  as  follows: 


Missile 

Targets 

Missile 

Targets 

1 

1,5 

6 

3,6 

2 

1,6 

7 

3,4 

3 

1,2 

8 

4,5 

4 

2,5 

9 

4,6 

5 

2,3 

We  will  represent  this  programming  of  missiles  by  a  Qualification  matrix 
Q: 


/iiioooooo\ 
[ooiiioooo 
000011100 
0  0  0  0  0  0  1  1  1 
IlOO  100010 
\o  1  0  0  0  1  0  0  l/. 


If  =  1,  this  means  that  missile  j  is  qualified,  i.e.  programmed, 

to  go  to  target  i.  Thus,  for  Q,  we  have  q^^  =  q^2  “  ^13  ”  so 
missiles  1,  2,  and  3  can  all  be  fired  at  target  1. 


Now  consider  the  situation  if  battle  actually  occurs.  It  may  well 
be  that  not  all  missiles  are  availablej  some  might  be  out  of  order, 
disable  by  enemy  action,  or  already  used.  Suppose  that  missiles  3,  4, 
and  8  have  been  eliminated.  Then  the  situation  is  represented  by  this 
qualification  matrix: 
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1  0  0  0  o\ 

0  10  0  0 
0  1110 
0  0  0  i  1 

0  0  0  0  0 

1  0  1  0  ly. 

The  first  problem  to  consider  is:  Find  an  optimal  assignment 
of  missiles  to  targets  given  a  qualification  matrix  such  as  We 

will  call  this  the  basic  assignment  problem.  In  this  basic  assignment 
problem  we  assume  that  all  missiles  are  to  be  launched j  we  do  not 
consider  the  possibility  of  keeping  some  missiles  in  reserve. 

Clearly,  the  first  thing  we  must  do  in  attacking  this  problem  is 
to  define  the  term  optimal.  We  propose  the  use  of  a  function  that  gives 
for  each  assignment  a  nmber  called  the  damage  for  that  assignment.  An 
optimal  assignment  is  then  one  for  which  the  damage  is  a  maximum.  A 
natural  idea  is  to  make  the  damage  equal  to  the  expected  value  of  the 
loss  which  the  enemy  will  suffer  if  the  missiles  are  launched  with 
the  given  assignment  of  targets.  However,  by  suitably  choosing  the 
damage  function,  many  different  optimization  criteria  may  be  represented. 
The  appendix  contains  several  examples  to  illustrate  this.  It  is  hoped 
that  these  examples  will  help  the  reader  gain  familiarity  with  the 
damage  function  concept.  Example  4»  which  shows  how  to  maximize  the 
probability  that  all  targets  are  destroyed,  is  especially  interesting. 
Precise  definitions  of  the  terms  "damage  function"  euid  "basic  assignment 
problem"  are  given  in  Section  1,  However,  no  discussion  of  the  solution 
of  the  problem  is  given.  For  that,  the  reader  must  turn  to  Reference  [l], 


•t 


which  he  should  certainly  read  along  with  this  document.  There 

an  efficient  algorithm  for  solving  the  basic  assignment  problem  is  described. 

The  second  problem  we  consider  is:  Evaluate  the  effectiveness  of 
the  original  qualification  matrix.  If  we  knew  that  all  missiles  would 
be  available  for  battle,  there  would  be  nothing  to  this  problem.  The 
measure  of  its  effectiveness  would  simply  be  the  damage  for  the  best 
assignment,  this  being  found  by  the  algorithm  mentioned  above.  However, 
we  do  not  expect  all  the  missiles  to  be  available  for  battle,  and  here 
a  random  element  enters  the  calculation.  It  seems  reasonable  to  assume 
that  each  missile  has  probability  p  of  being  available.  We  then  ask 
what  is  the  expected  value  of  the  damage  resulting  fl*om  the  best  possible 
assignment  of  the  available  missiles.  Except  in  some  simple  cases,  it 
seems  in^ossible  to  compute  this  expected  damage  by  analytic  methods.  Instead, 
we  turn  to  the  Monte  Carlo  method.  The  problem  could  be  attacked 
by  straightforward  simulation,  but  a  better  approach  is  possible.  This 
is  given  in  Section  2.  In  addition  to  asking  for  the  expected  damage, 
one  might  want  to  know  the  probability  that  the  damage  would  be  at  least 
a  certain  value.  For  example,  one  might  ask  for  the  probability  that 
all  targets  would  be  destroyed.  This  may  also  be  calculated  by  the 
Monte  Carlo  method.  The  calculation  is  described  in  Section  3. 

1,  Damage  functions  and  the  basic  assignment  problem 

In  this  section,  we  give  a  precise  statement  of  the  basic  assignment 
problem.  The  terms  defined  here  will  be  used  freely  in  the  following 
sections.  No  description  of  the  algorithm  for  solving  the  problem  is 
given  here.  The  reader  must  obtain  all  information  about  that  from 


r 


5 

Reference  [l]|  and  it  is  strongly  reconunended  that  he  do  that  before 
going  on  to  Section  2.  Examples  illustrating  the  use  of  the  damage 
function  will  be  found  in  the  appendix.  The  reader  might  also  wish 
to  look  at  some  of  these  before  continuing  to  the  discussion  of  Monte  Carlo 
methods  in  Sections  2  and  3. 

We  begin  with  some  definitions,  A  qualification  matrix  Q  is  a 
zero-one  matrix,  i.e.  a  matrix  in  which  every  element  is  either  zero  or 
une,  which  has  at  least  one  non- zero  element  in  each  column,  A  quali¬ 
fication  matrix  is  assumed  to  represent  the  programming  of  a  set  of 
missiles  as  described  in  the  introduction.  The  restriction  that  every 
column  has  a  non-zero  element  means  simply  that  every  missile  is  qualified 
for  some  target. 

An  assignment  matrix  A  belonging  i>o  a  qualification  matrix  Q 
is  a  zero-one  matrix  with  the  following  properties. 

(1.1)  If  q. ,  =  0,  then  a.,  =0. 

ij  ’  ij 

(1.2)  For  each  j  there  exists  exactly  one  1  such  that  =  1, 

An  assignment  matrix  is  supposed  to  give  the  targets  at  which  the 
available  missiles  will  be  launched.  Condition  (l.l)  means  that  a  missile 
can  be  launched  only  at  a  target  for  vdiich  it  is  qualified.  Condition 

(1.2)  says  first  that  a  missile  can  be  launched  at  only  one  target,  i.e. 
it  cannot  go  in  two  directions  at  once,  and  second  that  it  must  be 
launched  at  some  target.  This  last  condition  we  mentioned  in  the  intro¬ 
duction;  to  repeat;  "in  the  basic  assignment  problem  all  missiles  are  to 
be  laxinched  at  once;  none  are  to  ^  kept  in  reserve." 
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For  an  assignment  matrix  A,  we  let  r^(A)  = 
is  the  number  of  missiles  assigned  to  target  i. 


a 


ij 


,  i.e.  T^ik) 


Let  Q  be  an  n  x  m  qualification  matrix.  A  damage  function  for 
Q  is  a  function  g(i»J)  defined  for  all  i  with  1  ^  i  ^  n  and  all 
j  with  0  <  j  ^  m  which  has  the  following  properties: 


(1.3)  g(i,0)  =  0  for  all  1. 

(1.4)  g(i>j)  is  each  i  an  Increasing  function  of  J. 

(1.5)  g(i»J)  is  for  each  1  a  concave  function  of  J. 

If  we  let  =  g(i,j)  -  g(i,J  -  D  for  j  ^  1,  Condition  1.5  is 

equivalent  to 

(1.6)  If  j  ^  k,  then  &(i,k)  ^  ?)(i,j). 

The  quantity  g(i,J)  way  be  interpreted  as  the  expected  loss  which  the 
enemy  will  suffer  if  J  missiles  are  launched  at  target  i.  If  we 
assume  that  the  enemy’s  total  loss  is  the  sum  of  the  losses  for  each 
target,  then  the  expected  loss  to  the  enemy  resulting  from  an  assignment 
according  to  the  assignment  matrix  A  will  be 

n 

(1.7)  d(A)  =  Z  g(l,r.(A)). 

1=1 

This  leads  us  to  the  following  formulation  of  the  basic  assignment 
problem:  Given  a  nimllfi cation  matrix  Q  and  damage  function  g  for 
Q,  find  an  assignment  matrix  A  which  maximizes  the  damage  d(A)  given 


by  (1.7). 
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An  efficient  algorithm  for  the  solution  of  this  problem  is  described 
in  Reference  [l].  In  that  paper,  the  problem  is  described  in  terns  of 
men  and  tasks,  but  translation  into  the  terms  of  the  missile  problem  is 
easy.  For  the  term  "man"  in  [l],  substitute  "missile",  for  "Job" 
substitute  "target",  for  "output"  substitute  "damage".  The  output  or 
damage  in  [l]  is  given  by  a  set  of  functions  f2^(J)>»»*iT^(J)»  The 
relation  to  the  damage  function  g(i,J)  is  simply  f^(j)  =  g(i»J)* 

Interpreting  g(i»j)  ss  an  expected  loss,  the  Conditions  (1.3) 
and  (1.4)  are  obvious.  The  Condition  (1.5)  is,  however,  a  significant 
restriction  on  the  class  of  damage  functions  which  we  allow.  It  says, 
roughly,  that  for  each  target  the  expected  additional  damage  done  by., 
sending  one  more  missile  to  that  target  decreases  as  the  number  of  missiles 
assigned  to  that  target  is  increased.  It  is  possible  to  consider  damage 
functions  without  the  restriction  (1.5),  but  this  condition  is  necessary 
for  use  of  the  algorithm  of  Reference  [l] . 

2.  Expected  damage  for  a  given  qualification  matrix 

In  both  this  section  and  the  next,  we  vdll  let  denote  a  random 
subset  from  the  set  of  integers  l,2,...,m.  The  number  m  will  be 
assumed  to  be  given  and  fixed  throughout  the  argument.  Each 
has  independently  probability  p  of  being  in  X^.  Thus  p  is  a 
parameter  in  the  probability  distribution  of  X^.  The  number  of  elements 

in  X  will  be  denoted  by  a(X  ).  The  probability  that  o(X  )  =  r 

p  P  P 

will  bo  denoted  by  q^,  which  is  just  a  binomial  probability,  i.e. 

(2.1)  Prob(o(Xp)  =  r)  =  q^  =  (”)p^(l  -  p)”^ 
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Note  that  the  value  of  is  a  function  of  p,  although  this  Is  not 
made  clear  by  the  notation. 


The  purpose  of  this  section  and  the  next  is  to  investigate  the 
following  situation.  A  damage  function  g(i,j)  and  an  n  x  m 
qualification  matrix  Q  are  given.  However,  for  the  conjectured  attack 
each  missile  has  only  probability  p  of  being  available.  Thus,  the 
qualification  matrix  of  the  available  missiles  is  a  random  variable, 
Q(Xp),  where  Q(Xp)  is  the  matrix  obtained  from  Q  by  dropping  out 
those  columns  which  are  not  in  X^.  The  damage  for  an  optimal  assign¬ 
ment  belonging  to  Q(Xp)  will  be  denoted  by  d(Q(Xp)).  We  are 

interested  in  the  properties  of  the  random  variable  d(Q(X  )). 

P 


In  this  section  we  study  the  problem  of  estimating  E[d(Q(Xp))]. 
For  convenience,  we  will  write  E(Q,p)  =  E[d(Q(Xp))].  The  expected 
damage  E(Q,p)  is  a  natural  measure  of  the  effectiveness  of  a  quali¬ 


fication  matrix.  If  the  value  of  p  is  chosen  realistically,  we  might 
well  choose  between  two  qualification  matrices  and  by 

comparing  E(Q^,p)  and  ECQ^,?). 


Here  is  a  simple  example  of  calculating  the  expected  damage 
E(Q,p).  Take  Q  to  be  the  matrix 


Suppose  the  damage  function  is 

5 

For  any  random  variable  Y,  ELY]  will  denote  the  expected  value  of  Y. 
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g(i,0)  =0  all  i 

g(i,j)  =1  1  <  j,  all  i. 

That  is,  the  damage  is  just  the  number  of  targets  covered  by  the 

assignment.  Now  it  is  easy  to  see  that  d(Q(X  ))  =  a(X  ).  The 

P  P 

probability  that  X^  contains  exactly  k  missiles  is 

3 

Therefore,  E(Q,p)  =  Z  k(?)p^(l  -  p)^”^.  For  p  =  l/3  this  is 

k=0 

2  2  3 

E(Q,p)  =  0  +  1  X  3  X  (i)  X  (|)  +  2  X  3  X  (i)  x  (|)  +  3  x  1  x  (i)  x  1 

=  1.0. 

For  p  =  2/3,  this  is 

2  2  3 

E(Q,p)  =  0  +  1  X  3  X  (|)  X  (i)  +  2  X  3  X  {|)  x  (i)  +  3  x  1  x  (|)  x  1 

=  2.0. 


Unfortunately,  an  exact  calculation  of  the  expected  damage  can 
be  made  only  in  very  simple  cases  such  as  that  given  in  the  example. 

For  large  qualification  matrices  and  more  complication  damage  functions, 
the  calculation  must  be  made  by  the  Monte  Carlo  method.  A  sin.y''^ 


approach  would  be  to  carry  out  a  direct  simulation  of  the  problem. 
One  could  generate  a  sequence  of  Independent  random  sets  Xp,...,X^ 
and  form  the  estimate  ^d(Q(X^)).  The  expected  value  of  1 

is  clearly  E(Q,p).  The  variance  of  ?  could  be  estimated  by 


N 

2  [d(Q(xb)]^  -  T^). 
i=l  P 
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Two  objections  can  be  raised  against  this  procedure.  First  of 

all,  it  yields  an  estimate  of  E(Q,P)  only  for  one  value  of  p.  To 

find  E(Q,P)  for  a  different  value  of  p,  the  entire  calculation  would 

have  to  be  done  again,  using  the  new  value  of  p  in  generating  the 

random  sets  X^.  Secondly,  this  procedure  takes  no  accoiint  of  the 

operation  of  the  algorithm  which  finds  the  optimal  assignment  for  Q(xJ) 

and  computes  d(Q(X^)).  As  well  as  finding  d(Q(X^)),  the  algorithm 

p  P 

obtains  as  intermediate  results  the  values  d(Q*^(X^))  for  j  =  l,...,a(Xp), 

where  Q'^(X^)  is  the  matrix  containing  only  the  first  j  columns  of 
P 

Q(X^).  The  direct  simulation  procedure  makes  no  use  of  the  values 
d(Q‘^(Xp))  for  j<o(X^). 


A  procedure  which  avoids  both  of  these  objections  has  proved  quite 
efficient  in  practice.  It  makes  use  of  a  random  permutation  d  of  the 
integers  l,2,...,ra,  i.e.  a  permutation  chosen  at  random  from  all  the 
ml  permutations  of  the  integers  1,2,., ,,m.  For  the  permutation  d, 
we  let  Q(d)  denote  the  qualification  matrix  obtain  from  Q  by 
permuting  the  colximns  of  Q  according  to  d.  Thus,  if  the  elements  of 
Q(d )  are  q*^,  we  have  q|j  =  ^i(j(j)»  where  d(j)  is  the  integer 
into  which  the  permutation  d  takes  j.  For  j  ^  m,  we  let  (d ) 
denote  the  n  x  j  qualification  matrix  consisting  of  the  first  j 
colimms  of  Q(d), 


To  estimate  E(Q,p)  proceed  as  follows;  Generate  a  sequence 
d  ^  ,d  2 , . .  •  »d  of  independent  random  permutations  of  the  integers 

1,2,.,, ,m • 

5  =  ^  Z  d(QJ(d.)). 

J  ^  i=l  ^ 


(2.2) 
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Then  for  any  probability  p,  an  unbiased  estimate  of  E(Q,p)  ^ 


m 


(2.3)  C?  -  ^  q.S., 
j=0  J  J 

where  qj  Is  given  by  (2.1). 


Note  that  the  parameter  p  enters  into  the  calculation  only 
through  the  q^  in  (2.3).  Thus,  the  number  of  calculations  involving 
the  parameter  p  is  independent  of  the  nimber  of  permutations  used. 
Moreover,  one  may  save  the  values  of  the  5^.  Then,  at  some  later  date, 
the  estimate  W  may  be  calculated  for  various  values  of  p  vd.thout 
generating  any  additional  random  permutations. 


The  natural  measure  of  the  error  in  a  Monte  Carlo  estimate  such 
as  (2.3)  is  the  standard  deviation  of  the  estimate.  This  is  defined 
as  follows.  The  variance  of  the  random  variable  W  is  V[W]  = 

E[  (^7  -  E[i?])^].  The  standard  deviation  of  R  is  .  The  usual 

procedure  in  a  Monte  Carlo  calculation  is  to  make  an  unbaised  estimate 
of  V[W]  and  then  take  the  square  root  of  this  as  an  estimate  of  the 
standard  deviation.  Therefore,  in  the  remainder  of  this  paper  we  will 
only  discuss  the  problem  of  estimating  the  variance. 


Unfortunately,  if  we  desire  to  estimate  the  variance  of  W,  seems  that 
intermediate  computations  Involving  p  must  be  made.  Let 


m 


(2.4)  ^i=  ^  qj,d(Q^(cJ^)), 

r=0 


and  let 

,  N 

(2.5)  ^  ^Z^d(Q'‘((5^))d(Q®(d^)). 

An  unbiased  estimate  of  V[Vl]  may  be  computed  by  using  either  of  the 
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followlnr  two  expressions 

(2.6)  z  =  n^n  -  ^  ^ 

,  m  « 

I*  I S  “O 

The  two  expressions  on  the  rifrht  are,  in  fact,  equal. 


W^] 


or 


The  expression  (2.7)  might  be  used  to  avoid  intermediate  computations 
with  the  parameter  p.  Unfortunately,  many  computations  vd.ll  be  needed 
to  find  the  G  ,  In  fact,  it  vn.ll  require  m(m  -  l)/2  multiplications 
and  additions  for  each  randan  permutation. 


The  remainder  of  this  section  will  be  devoted  to  proving  that  the 
estimates  and  2  do  have  the  desired  properties.  Note  that  we  can 
vn*ite 

,  N  m  ,  N 

{2.8)  0  =  i  2  [  2  q  d(Q''(o,))]  =5  2  W  . 

^  i=l  r=0  ^ 

Thus,  is  the  average  of  the  identically  distributed  random  variables 
Wi,W2,...,Wj^.  Thus  we  have 

m 

(2.9)  E[B]  =  E[Wj  =  2  q  E(<i(Q''{a))] 

^  r=0  ^ 

where  6  is  a  random  permutation.  Moreover,  by  familiar  argument  from 
the  theory  of  statistics,  we  know  that  an  unbiased  estimate  of  V[W]  is 
just  Z  as  given  by  (2.6).  To  verify  that  E[il]  is  equal  to 
E[d(Q(Xp))],  we  need  only  note  that 

(2.10)  E[d(Q^(d))]  -  E[d(Q(Xp))la(Xp)  =  r] . 
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Here,  the  term  on  the  right  denotes  the  conditional  expectation  of  the 

random  variable  d(Q(Xp))  given  that  =  r,  i.e,  that  the  number 

of  elements  in  X  is  equal  to  r.  The  expression  on  the  right  of 

P 

(2.10)  is  independent  of  p,  because  all  sets  with  the  same  number  of 

elements  have  equal  probability.  Now  by  definition  is  equal  to 

Prob[a(Xp)  =  r] .  Thus  (2,9)  becomes 

m 

(2.11)  E[^J]  =  Z  E[d(Q(X  ))Ia(X  )  =  r]Prob[a(X  )  =  r]  =  E[d(Q(X  ))],  q,e,d, 

r=0  ^  P  P  P 

The  only  remaining  statement  to  prove  is  that  (2,6)  and  (2,7)  are 
equivalent,  i,e,  that 

,  N  p  m 

(2.12)  ^  2  W,  =  2  q  q  G  . 

N  .  T  i  ^  rs 

1=1  r,s=0 

This  is  obtained  by  squaring  the  expression  (2.4) •  We  get 
,  N  ,  .  N  ra 

i=l  i=l  r,s=0 

m  T  N 

=  Z  q  q^[i  Z  d((3''f<j^))d(5®{Oj))] 

r,3=0  ^  i=l 


r,s=0 

3.  Probability  that  the  damage  exceeds  a  given  value 

The  notation  of  the  first  two  paragraphs  of  Section  2  applies  also 

in  this  section.  We  wish  to  study  further  the  random  variable  d(Q(X  )), 

P 

The  expected  value  E(Q,p),  discussed  in  Section  2,  is  certainly  the  most 
useful  single  number  associated  with  the  random  variable.  It  is  possible, 
however,  to  get  more  detailed  information  about  the  distribution  of 


d(Q(Xp)).  For  any  number  c,  we  can  estimate  Prob[d(Q(Xp) )  >  c] ,  i.e. 


u 


the  p’’obablllty  that  d(Q(Xp))  will  be  at  least  c.  Knowing 

Prob[d(Q(X  ))  2  c]  for  several  values  of  c  would  give  us  a  good  idea 
P 

of  the  distribution  of  the  random  variable  d(Q(Xp)).  In  some  cases, 
knowing  Prob[d(Q(Xp) )  ^  c]  for  a  single  value  of  c  would  be  very 
useful.  For  example,  if  the  damage  is  Just  the  number  of  targets  covered 
by  the  assignment  and  c  is  the  total  number  of  targets,  then 
Prob[d(Q(Xp) )  ^  c]  is  the  probability  of  covering  all  targets. 

With  notation  the  same  as  in  Section  2,  we  can  describe  the  Monte 
Carlo  procedure  to  estimate  Prob[d(Q(Xp))  ^  c]  as  follows: 

Let  the  damaf^e  function  g(i,j),  the  qualification  matrix  Q,  and 
the  number  c  be  given.  Generate  a  sequence  of  independent  random 
permutations  For  1  ^  i  ^  N  and  0  <  r  <  m,  set 

=  1  If  d(Q’'(d^))  2  0, 

(3.1) 

F^(aj)  =  0  if  d{Q''(a^))  <  c. 
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Note  that  for  fixed  i,  d(Q^(c3^))  ^  •••  d(Q''^(c3^))»  'fherefore,  for 

each  permutation  o^,  there  is  a  critical  number  n(c)^)  such  that 

FrC^i)  =  0  r  <  n(o^) 

(3.6) 

Fr(c^i)  =  1  r  ^  n(d^). 

To  determine  F^(d^)  for  all  r,  it  is  only  necessary  to  find  n(c5^). 

In  order  to  verify  that  (3.3)  is  indeed  an  unbiased  estimate  of 
Prob[d(Q(Xp))  >  c],  note  that 

m 

(3.7)  E[?]  =  2  q^E[F  (d)],  where 

r=0  ^  ^ 

d  is  a  random  permutation.  Further,  we  have 

(3.8)  E[F  (d)]  =  Prob[d(Q(X  ))  ^  c|a(X  )  =  r] 

i  Jr  H 

and 


(3.9) 


Prob(Q(X  )  =  r)  =  q^, 


(Recall  that  ci(Xp)  is  the  number  of  missiles  in  X^j  and  that  given 
that  ci(Xp)  is  equal  to  r,  X^  is  simply  a  random  subset  of  r 
elements,  all  such  subsets  being  equally  probable.)  Combining  (3.7)  - 

(3.9) ,  we  have 

m 

(3.10)  E[Y]  =  Z  Prob(a(X^)  =  r )Prob[d(Q(X^) )  ^  cla(X  )  =  r] 

r=0  ‘  P  P 


=  Prob[d(Q(Xp) )  ^  c],  q.e.d. 


16 


To  see  that  (3.4)  an  unbiased  estimate  of  V[Y] i  let 


r=0 

Then  the  are  independent  identically  distributed  random  variables, 

and  ?  -  “2  Y^,  Therefore,  just  as  in  expression  (2,6),  we  have  that  an 
unbiased  estimate  of  V[?]  is 

'3-11)  z  =  - ’"l  • 

2 

To  prove  that  this  is  equivalent  to  (3*4) >  we  expand  Y^ 

(3.12)  Yf  =  (  Z  q  F  (d  ))'"  =  Z  [qV(dj)  +  22  ^  * 

r=0  r=0  s=r+l 

Now  F  (d.)  is  either  zero  or  one,  and  for  r  <  s,  F  (d.)  ^  F  (d.), 
r  i  '  r  i  s  i 

Therefore,  (3.10)  becomes 

P  m  p  m  m 

(3.13)  Y^=  "['0^2  V^Vh>- 

r=0  s=r+l  r=0 

Thus,  we  have 


1  ^  P  ^ 

^  2  Yf  =  2  t  F  , 
N  1  r  r' 


i=l  ^  r-O 


and  this  proves  that  (3.11)  is  equivalent  to  (3.4) >  q.e.d. 
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Appendix 

In  this  appendix  we  present  four  examples  to  illustrate  the  use 
of  the  damape  function  idea.  The  first  example  pives  the  natural 
construction  of  a  damage  function  based  on  the  enemy's  expected  loss. 

The  next  two  examples  show  how  the  damage  function  may  be  used  to  define 
an  assignment  problem  based  on  the  idea  of  target  priorities.  This 
includes  the  simple  special  case  of  assigning  one  missile  to  each  target. 
The  last  example  deals  with  the  problem  of  maximizing  the  probability 
that  all  targets  be  destroyed.  This  is  easily  done  using  a  suitable 
damage  function  and  the  algorithm  of  Reference  [l],  yet  it  does  not 
appear  on  the  surface  to  be  a  problem  accessible  to  the  damage  function 
approach. 

Example  1.  Damage  function  by  analysis  of  target  values  and  kill 
probabilities 

There  are  six  targets,  which  have  the  following  values  to  the 
enemy: 


Targets  Value  each  target 

1,  2,  3  20 

4,  5,  6  10. 

The  value  is  not  assixmed  to  be  in  any  concrete  measure  such  as  men  or 
dollars,  although  it  may  be;  it  is  simply  a  measure  of  the  relative 
importance  of  the  target.  Thus  target  2,  with  a  value  of  20,  is  twice 
as  important  as  target  5,  which  has  a  value  of  only  10.  From  this,  we 
conclude  that  an  action  which  would  destroy  target  2  with  probability 
1/2  is  equivalent  to  one  which  would  destroy  target  5  with  certainty. 
In  each  case,  the  expected  loss  to  the  enemy  is  10, 
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To  conotruct  the  damage  function,  it  is  necessary  to  have  the  kill 
probability  for  each  target,  i.e*  the  probability  that  a  single  missile 
will  destroy  the  target.  The  next  table  gives  these  kill  probabilities. 

Targets  Kill  probability 

1,6  *2 

2,  5  .5 

3,  4 

We  assume  that  if  several  missiles  are  fired  at  a  target  their 
effects  are  independent.  Each  missile  either  destroys  the  target  or 
leaves  it  undamaged.  Then,  if  p^  is  the  probability  that  a  single 
missile  will  destroy  target  i,  the  probability  that  a  salvo  of  j 
mis  lies  will  destroy  target  1  is 

1  -  (1  - 

Thus,  if  the  value  of  target  i  is  V^,  the  expected  damage  to  the  enemy 
from  a  salvo  of  J  missiles  at  target  i  is 

g(i,J)  =  Vj[l  -  (1  - 

Note  that 

=  g(l.J)  -  gd.J  -  1)  =  (1  -  Pi)^‘^(VjP^). 

Thus  f)(i,j)  decreases  with  j;  g(i,j)  is  concave  in  j. 

The  next  table  gives  the  value  of  g(i,j)  for  j  ^  6.  Values 
for  J  greater  than  6  could  be  easily  computed  from  the  information 
given  above. 
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Table  1 


i/j 

0 

1 

2 

3 

4. 

5 

6 

1 

0 

4 

7.2 

9.76 

11.81 

13.46 

14.76 

2 

0 

10 

15 

17.5 

18.75 

19.38 

19.69 

3 

0 

16 

19.2 

19.84 

19.99 

20,00 

20.00 

4 

0 

2 

3.6 

4.88 

5.90 

6.72 

7.38 

5 

0 

5 

7.5 

8.75 

9.38 

9.69 

9.84 

6 

0 

8 

9.6 

9.92 

9.99 

10.00 

10.00 

The  value  of  the  damage  for  any  assignment  matrix  may  be  easily 
computed  from  this  table.  Suppose,  for  example,  we  have  nine  missiles 
with  this  qualification  matrix  Q 


A 

1 

1 

0 

0 

0 

0 

0 

A 

0 

1 

1 

1 

0 

0 

0 

0  ' 

Q  = 

0 

0 

0 

0 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

0 

0 

1 

0 

0 

0 

1 

°l 

\o 

1 

0 

0 

0 

1 

0 

0 

1 

• 

The  following  assignment  matrix  A  belonging  to  Q  assigns  two  missiles 
to  each  of  the  targets  with  value  20  (numbers  1,  2,  and  3)  and  cne  to 
each  of  the  other  targets, 

A  10000000 
001100000 
,  000011000 
000000100 
000000010 
yo  00000001 

For  i  =  1,  2,  and  3,  r^(A)  =  2;  and  for  i  =  4,  5,  and  6,  r^(A)  =  1, 

Thus  the  damage  d(A)  is  given  by 
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6  3  6 

d(A)  =  Z  f;(i,r,(A))  =  Z  g(i,2)  +  Z  c(i»l) 

i=l  i=l  1=4 

=  7.2  +  15.0  +  19.2  +  2.0  +  5.0  +  8.0  =  56.4. 

Is  this  an  optimal  assignment  for  Q?  This  is  not  easy  to  determine 
using  Table  1.  Much  mo.  e  convenient  for  this  purpose  is  a  table  of  the 
difference  The  difference  table  corresponding  to  Table  1  is: 


Table  2 


i/J 

1  2 

3 

4 

5 

6 

1 

4.00  3.20 

2.56 

2.05 

1.65 

1.30 

2 

10.00  5.00 

2.50 

1.25 

.63 

.31 

3 

16.00  3.20 

.64 

.15 

.04 

.00 

4 

• 

2.00  1.60 

1.28 

1.02 

.82 

.66 

5 

5.00  2.50 

1.25 

.63 

.31 

.15 

6 

8.00  1.60 

.32 

.07 

.01 

.00 

The 

next 

table  reproduces  the 

first  3  coliimns  of  Table  2, 

In  additii 

certain 

entries  are  starred. 

The  number 

of  starred 

entries 

in  a  row 

is 

the  number  of  missiles  assigned  to  that  target  by  the  matrix  A, 

Table  3 

i/j 

1 

2 

3 

1 

4.00»‘ 

3.20» 

2.56 

2 

10.00^^ 

5.00* 

2.50 

3 

16.00» 

3.20* 

.64 

4 

2.00^ 

1.60 

1,28 

5 

5.00^^ 

2.50 

1.25 

6 


8.00^'‘ 


1.60 


.32 
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ITie  damage  for  the  aaaif'nment  A  is  Just  the  num  of  the  starred  entries 
in  Table  3.  Now  look  at  the  two  underlined  entries  *,  (4,1)  "  2,00  and 
^(1,3)  =  2,56.  The  value  of  the  assiprinent  matrix  would  clearly  be  hi^’her 
if  t(l,3)  was  starred  instead  of  ^(4,l).  An  assignment  which  effects 
this  improvement  is: 


1 

0 

0 

0 

0 

0 


1 

0 

0 

0 

0 

0 


0 

1 

0 

0 

0 

0 


0  0  0 
10  0 
oil 
0  0  0 
0  0  0 
0  0  0 


0 

0 

0 

0 

1 

0 


The  damage  for  B  is  56.96,  which  is  .56  greater  than  that  for  A,  i.e. 
just  the  difference  between  t(l,3)  and  ^>(3,1).  The  assignment  matrix 
B  is  in  fact  an  optimal  assignment  matrix  belonging  to  Q.  Moreover, 
no  assignment  matrix  for  nine  missiles  can  improve  on  B  whatever  the 
qualification  matrix  may  be.  To  see  this,  simply  note  that  in  Table  3 
with  J)(l,3)  starred  instead  of  ^(4,l)  the  nine  largest  values  are 
starred.  Now  the  damage  for  any  assigrLment  of  nine  missiles  will  be  a 
sum  of  nine  terms  taken  from  Table  2.  Since  the  dairiage  for  B  is  the  rum 
of  the  nine  largest  terms,  it  is  the  maximum  possible. 


An  Interesting  point  to  note  about  Table  2  is  that  ^(3,3)  is 
less  than  t(4,3)  even  though  the  value  of  target  4  is  only  half  that 
of  target  3.  This  occurs  because  target  3  is  relatively  soft  and  will 
almost  certainly  be  destroyed  by  a  salvo  of  two  missiles.  In  fact,  the 
probability  of  the  missiles  destroying  target  3  is  .96,  so  there  is  very 
little  to  gain  by  adding  a  third  missile  to  this  target. 


hxar.ple  . 


P? 


On 


*'  la 


Inslle  to  *‘aGh  tar/^et  ^Ith  a  i;et  ol'  target  priorities 


ri'.f'  probleia  ol‘  assipnln,’  onf^  rniiisile  to  each  target  i  r>  n  special 
ca.'e  ol’  the  basic  assipnr.ent  problem,  Suppose  there  are  six  missiles 
anJ  six  tar-'ots  with  the  ‘‘ollowinw  quali  I'i  cat:  on  matrix 


Q  = 


/:  0  0  0  c  0^ 

l  0  0  0  0  0 

0  110  10 
0  10  10  1 

0  0  0  1  1  0 

\  0  0  1  0  0  1  / 

\l  1  1  1  1  1/ 


Tiiere  is  an  extra  row  in  Q,  This  has  been  introduced  as  a  dummy 
target,  for  which  all  missiles  are  qualified,  to  take  care  of  missiles 
not  otherwise  assigned.  It  is  necessary,  because  at  most  one  missile 
is  to  be  assigned  to  each  real  target. 


Suppose  that  it  is  desired  to  find  an  assignment  which  assigns  at 
most  one  missile  to  each  of  the  targets  1-6  and  covers  as  many  of  these 
targets  as  possible.  If  the  following  damage  function  is  used,  any 
optimal  assifumient  for  the  basic  assif’^nment  problem,  will  have  the  desired 
properties. 

g(i,0)  ^0  all  i 

(1)  g(i,j)  ^1  1  ^  i  <'  6,  1  <  j 

f(7,j)  -  j/P  1  J  . 

Note  that  ^(i,j)  -  0  for  1  6  and  j  ^  P,  and  tCTjj)  =  l/P  for  all 


j.  Therefore,  because  all  missiles  are  qualified  for  target  7,  an  optimal 
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a33if»nment  will  a33i(’:n  at  most  one  missile  to  each  of  the  targets  1-6. 

If  n  missiles  are  assigned  to  n  of  the  targets  1-6  and  the  remainder 
to  target  7,  the  dummy  target,  the  damage  vd.ll  be  n  +  l/2(6  -  n). 
Clearly  this  vfill  be  a  maximum  v;hen  n  is  a  maximum.  Thus  any  optimal 
asrdgnment  has  the  desired  properties. 


Here  are  tvro  optimal  assignment  matrices  for  this  problem 


(2) 


A  ^ 


t  ° 

/  0  0 

0  1 
0  0 
0  0 
0  0 
0  0 


0  0  0 
0  0  0 
0  0  0 
0  10 
0  0  1 
0  0  0 
10  0 


0  0 
0  0 
1  0 
0  0 
0  0 
0  0 
0  1 


0  0 
0  0 
0  0 
1  0 
0  1 
0  0 
0  0 


The  assignment  A  fails  to  cover  target  2  v/hile  B  fails  to  cover 
target  1.  It  is  not  possible  for  an  assignment  belonging  to  Q  to 
cover  both  of  these  targets. 


It  is  possible  to  constinict  a  damage  function  for  which  any  optimal 
matrix  will  have  the  properties  given  above,  and  in  addition,  the 
assignment  matrix  A  vdll  be  preferred  over  B,  That  is  to  say.,  a  set 
of  priorities  can  be  introduced  for  the  targets  so  that  target  1  is 
more  important  than  target  2.  To  give  an  exact  statement  of  a  priority 
scheme,  it  is  necessary  to  give  an  ordering  for  the  assignment  matrices. 
An  optimal  matrix  belonging  to  0  is  then  one  which  is  greater  than  all 
other  matrices  belonging  to  Q. 

Suppose  that  target  1  is  most  important,  target  2  is  next  most 
important,  and  so  on  to  target  6,  the  least  important  target.  An 
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asnif^nment  problem  for  thio  net  of  targets  and  priorities  might  be 
stated  in  this  way: 


Priority  Assignment  Problem 

Consider  only  m*atrices  belonging  to  Q  for  which 
r^(A)  <  1  for  1  <  i  <  6.  For  two  such  matrices  A  and  B, 
we  say  that  A  is  greater  than  or  equal  to  B  and  uTite 
B  <  A  if  one  of  the  following  conditions  holds: 

6  6 

(i)  Z  r.  (B)  <  4  r,  (A), 

i=l  i=l  ^ 

6  6 

(ii)  Z  r,  (B)  =  Z  r,  (A)  and 

i=l  ^  i=l  ^ 

for  some  k  <  6,  r^(B)  <  while  rj^(A)  =  r^(B) 

for  all  1  <  k, 

(iii)  r^(B)  =  r^(A)  for  all  i  <  6. 

The  problem  is  to  find  an  assignment  matrix  A  belonging  to 
Q  and  with  r^(A)  <  1  for  i  <  6  such  that  any  other  matrix 
B  with  these  properties  satisfies  B  <  A. 


Consider  the  assignment  matrices  A  and  B  of  (2).  The 

6 

assignment  matrix  A  is  optimal.  We  have,  in  fact,  Z  r  (A)  =  5, 

i=l  ^ 

r^(A)  =  1,  r^{A)  =  0,  r^(A)  =  1  for  all  i  with  3  <  i  <  6.  It  was 

pointed  out  above  that  no  matrix  belonging  to  Q  could  assign  to  both 

6 

targets  1  and  2,  Hence,  the  maximum  possible  value  of  Z  r.  (A)  is 

i=l  ^ 

5;  and  if  r^(A)  =  1,  r2(A)  must  be  zero.  The  matrix  A  is  greater 
than  B,  because  r^(A)  -  1  while  r^(B)  =  0. 
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A  damage  function  for  which  any  optimal  matrix  will  be  a  solution 
to  the  target  priority  problem  is  easily  constructed  by  modifying  the 
damage  function  (l).  For  the  function  of  (l),  the  quantities 
f)  (i  ,1) ,  1  <  1  <  6  are  all  equal  to  one.  It  is  only  necessary  to  make 
these  increases  in  order  of  target  priority,  provided  that  we  keep 
?)(i,l)  >  h(7,l)  for  1  <  i  <  6.  For  example,  we  can  use 

g(i,0)  =0  all  i 

(4)  e(l,J)  =  7-1  1<1<6,  1<J 

s(7,j)=j/2  1  <  J. 

The  fact  that  a  damage  function  of  this  type  solves  the  priority  assign¬ 
ment  problem  is  proved  in  [l]. 

Example  3>  A  more  complicated  set  of  priorities 

Consider  the  damage  function  g(i,j)  whose  differences  f)(i,j) 
are  given  by  this  matrix 


/6.0 

6.0 

0.0\ 

f  5.0 

5.0 

0.0  \ 

4.0 

2.0 

0.0 

3.0 

0.0 

0.0 

1 

0.0 

0.0  ) 

0.0 

o.oy 

Write  the  targets  in  a  list  entering  target  i  in  the  list  one 
time  for  each  non-zero  b(i,J),  Start  with  the  largest  ?)(i,j)  and 
work  downward.  The  result  is  this  list 


(6) 


2,  2,  3,  4>  3,  5f  6, 
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Suppose  that  n  missiles  are  available,  eacn  qualified  for  all  targets. 
An  optimal  assignment  for  the  damage  function  (5)  will  assign  missiles 
to  the  first  n  targets  in  the  list  (6),  each  target  getting  one  missile 
for  each  of  its  occurences  in  the  first  n  elements  of  the  list.  For 
example,  if  five  missiles  are  available,  they  will  be  assigned  to  targets 

1,  1,  2,  2,  and  3,  i.e.  two  missiles  to  target  1,  two  missiles  to  target 

2,  and  one  missile  to  target  3,  With  seven  missiles  available,  targets 
1,  2,  and  3  would  each  receive  two  missiles,  and  one  missile  would  be 
assigned  to  target  4, 

When  the  missiles  are  not  qualified  for  every  target,  the  assignment 
will  in  a  certain  sense  come  as  close  as  possible  to  assigning  missiles 
to  the  first  n  targets  in  the  list.  To  be  more  precise,  an  optimal 
assignment  A  will  have  the  following  properties: 

(i)  If  possible  A  will  assign  two  missiles  to  target  1,  If 
that  is  not  possible  and  one  missile  can  be  assigned  to  target  1,  then 
it  will  be, 

(ii)  In  addition  to  satisfying  (i),  A  will  assign,  if  possible,  two 
missiles  to  target  2,  If  that  is  not  possible  but  one  missile  can  be 
assigned  to  target  2,  then  it  will  be. 

(iii)  In  addition  to  satisfying  (i)  and  (ii),  A  will  if  possible 
assign  a  missile  to  target  3. 

(iv)  In  addition  to  satisfying  (i)  -  (iii),  A  will  if  possible 
assign  a  missile  to  target  4. 

(v)  In  addition  to  satisfying  (i)  -  (iv),  A  will  if  possible  assign 
a  second  missile  to  target  3, 

(vi)  and  so  on. 
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The  relationship  between  this  set  of  conditions  and  the  list  of 
tarf^ets  (6)  should  be  clear  to  the  reader.  We  call  the  list  (6)  a 
list  of  priorities  and  say  that  an  assignment  matrix  which  is  optimal 
for  the  damage  function  (5)  is  optimal  with  respect  to  this  list  of 
priorities.  The  list  of  priorities  is  not  merely  an  ordering  of  the 
targets.  It  specifies  how  many  missiles  should  be  assigned  to  each 
target  and  in  what  order  the  possibly  conflicting  requirements  should 
be  met.  For  example,  with  the  list  (6),  one  missile  will  be  assigned 
to  target  3  rather  than  target  4.  However,  assigning  a  missile  to 
target  4  is  preferred  over  assigning  the  second  missile  to  target  3* 
Moreover,  no  preference  is  made  between  target  5  and  6.  Let  us  see  how 
all  this  will  work  out  with  a  qualification  matrix  Q, 

/l  1111100  o\ 

000011111 

1^^110000100 

001000010 

000100000 

\ooooooooi/. 


An  optimal  assignment  matrix  belonging  to  Q  with  the  damage 
function  (5)  is 


A 

0 

0 

0 

1 

1 

0 

0 

A 

1  0 

0 

0 

0 

0 

0 

1 

1 

0^ 

1 

1 

0 

0 

0 

0 

0 

0 

0 

A  = 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

\o 

0 

0 

0 

0 

0 

0 

0 

1 

• 

This  achieves  the  maximum  possible  damage  of  34. 
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If  missiles  1 
matrix  is: 

A  ' 

0  0 


,  ", 


and  7  of  0  are  not  available,  the  qualification 


«2  = 


0  0 
1 
0 


1 
1 
0 

0  0 
1  0 


V 


0  0  0  0 


An  optimal  assif^nment  matrix  for  this  is 


/ 


A.  = 


1  1  0  0  0  0 
0  0  110  0 
0  0  0  0  0  0 
0  0  0  0  1  0 
0  0  0  0  0  0 
0  0  0  0  0  1 


\ 


The  damage  for  this  assignment  is  26,5. 


If  missiles  3,  4,  8,  and  9  are  not  available  from  Q,  the  qualification 
matrix  is 


Q.  = 


An  optimal  assignment  matrix  belonging  to  0^  is 
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M 


3 


1 

0 

0 

0 

0 

0 


The  damaf^e  for 


0  0 
1  1 
0  0 
0  0 
0  0 
0  0 

is  26. 


Now  as  is  indicated  in  [l],  the  magnitudes  of  the  ^(i,j)  are  not 
important  in  determining  the  list  of  priorities.  Only  the  ordering  of 
the  is  important.  Thus,  the  damage  function  g^,(i,j)  with 

different  ^^(i,j)  given  belo'*,  results  in  the  same  optimal  assignments 
as  g(i,j)  above. 


20  o\ 
17  0  ] 
2  0 
0  0 
0  0 
0  0 


However,  for  the  assignment  matrices  and  A^,  the  damages  using 

g^(i,.')  are  79  and  90  respectively.  Thus,  using  g„(i,.')  it  seems 
better  to  have  missiles  1,  2,  5,  6,  and  7  available  (qualification 
matrix  Q^)  than  to  have  missiles  3>  4>  5,  6,  8,  and  9  available 
(matrix  ^2^*  other  hand,  the  reverse  is  true  if  we  use  the 

damage  function  g(i,j). 


The  point  of  all  this  is  that  it  is  not  really  sufficient  to 
specify  a  list  of  target  priorities.  If  we  are  to  be  able  to  compare 
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the  efrectiveneas  of  different  qualification  matrices,  we  must  have 
more  indication  of  the  relative  values  of  the  various  targets. 
Constructing  the  damage  function  by  a  procedure  like  that  of  Example  1 
will  give  such  information. 

Example  4.  Maximizing  the  probability  of  destroying  all  targets 

Assume  that  for  each  target  1  the  kill  probability  p^  Is  known. 

Then  the  probability  that  target  1  will  be  destroyed  if  missiles  are 

r^(A) 

launched  according  to  an  assignment  matrix  A  Is  l-(l-p^)  • 

Therefore,  the  probability  that  all  targets  will  be  destroyed  is 

n  r  (A) 

(7)  P.  =  n  (1  -  (1  -  p  )  ). 

^  i^l  ^ 

Given  a  qualification  matrix  0,  we  ask  for  an  assignment  matrix  A 
belonging  to  Q  for  which  P(A)  is  maximum. 

At  first  it  might  appear  that  this  problem  cannot  be  formulated 
using  a  damage  function  as  described  in  Section  1.  However,  note  that 
P(A)  is  a  maximum  if  and  only  if  log(P(A))  is  a  maximum,  and 

n  r  (A) 

(8)  log  P(A)  =  2  log(l  -  (1  -  p  )  ^  ). 

1=1  ^ 

Now  this  looks  very  much  like  a  damage  function  expression.  Let  us 
define  g(i,j)  as  follows; 

g(i,0)  =0  all  i 

g(i,j)  =  c  +  logd  -  (1  -  p^)*^')  for  j  ^  0, 
where  c  is  chosen  so  that  for  all  i  and  k 
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(9)  c(i>l)  ^  -  c(k,l). 

We  claim  that  ic  a  damage  fimction,  i.e,  it  satisfies 

(1.3)  -  (1.5).  Clearly  g(l,0)  =  0.  To  verify  that  g(i,j)  is  concave 
in  j,  let  f(t)  =  (1  -  (1  -  P^)^)  Tor  t  >  0  and  h(t)  =  logf(t). 

In  Example  1  it  was  proved  that  f  is  concave,  i.e.  that  f”(t)  ^  0. 
Clearly  f’(t)  ^  0  and  f(t)  ^  0.  Therefore,  for  h(t),  we  have 


(10) 

and 

(11) 


h»(t)  =  f»(t)/f(t)  ^  0 


h"(t) 


<  0. 


From  (lO),  we  conclude  that  for  all  i  and  j  ^  1»  g(i>j)  is  increasing 

in  j,  and  from  (11),  we  conclude  that  for  all  i  and  J  >  1,  g(iJ) 

is  concave  in  j.  The  choice  of  c  so  that  (9)  holds  insures  that 

g(i,j)  is  concave  and  increasing  in  j  for  all  j  ^  0. 


Now  suppose  that  it  is  possible  to  find  an  assignment  matrix  A* 
belonging  to  Q  such  that  for  all  i,  r^{k^)  ^  0.  Then  we  claim  that 
log  P(A)  is  a  maximiim  if  and  only  if 

n 

d(A)  =  Z  g(i,r(A)) 
i=l 

is  a  maximum.  To  see  this,  simply  note  that  because  of  Condition  (9) 
d(A)  cannot  be  a  maximum  if  r^(A)  =  0  for  some  i.  But  with 
r.(A)  f  0  for  all  i,  we  have 

d(A)  =  log  (P(A))  +  nc. 

Now  if  no  matrix  A  belonging  to  Q  exists  for  which  r^(A)  ^  0 
for  all  i,  then  P(A)  =  0  for  all  matrices  A  belonging  to  Q.  Hence 


wo  can  conclude  that  f'or  an  anni,'^nment  matrix  A  bolonf^ln,'^  to  0, 
F(A)  in  a  m'ucinum  it  the  ioj^ai^e  d(A)  in  a  maximum# 
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